Classical aspects of quantum fields far from equilibrium 
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We consider the time evolution of nonequilibrium quantum scalar fields in the 0{N) model, using 
the next-to-leading order 1/A'' expansion of the 2PI effective action. A comparison with exact 
numerical simulations in 1 -I- 1 dimensions in the classical limit shows that the 1/A expansion 
gives quantitatively precise results already for moderate values of N. For sufficiently high initial 
occupation numbers the time evolution of quantum fields is shown to be accurately described by 
classical physics. Eventually the correspondence breaks down due to the difference between classical 
and quantum thermal equilibrium. 
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In recent years we have witnessed an enormous increase 
of interest in the dynamics of quantum fields out of equi- 
librium. Much progress has been achieved for systems 
close to thermal equilibrium or with effective descriptions 
based on a separation of scales in the weak coupling limit 
P[. Current and upcoming relativistic heavy-ion colli- 
sion experiments provide an important motivation to find 
controlled nonperturbative approximation schemes which 
yield a quantitative description of far-from-equilibrium 
phenomena from first principles. 

Practicable nonperturbative approximations may be 
based on the two-particle irreducible (2P/) generating 
functional for Green's functions [g[. Recently, a system- 
atic 1/N expansion of the 2PI effective action has been 
proposed for a scalar 0{N) symmetric quantum field the- 
ory H . This nonperturbative approach extends previous 
successful descriptions of the large-time behavior of quan- 
tum fields |3M| , which employ the loop expansion of the 
2PI effective action relevant at weak couplings |g,^ . At 
next-to-leading order (NLO) the 1/N expansion of the 
2PI effective action has been solved for the quantum 
theory in 1 -I- 1 dimensions |^]. The approach overcomes 
the problem of a secular time evolution, which is encoun- 
tered in the standard 1/A^ expansion of the IPI effective 
action beyond leading order M. 

Our purpose in this Letter is to establish that the 1/A^ 
expansion at NLO gives quantitatively precise results al- 
ready for moderate values of N. We therefore have a 
small nonperturbative expansion parameter at hand and 
a controlled description of far-from-equilibrium dynamics 
becomes possible. As an application we compare quan- 
tum and classical evolution and demonstrate that the 
nonequilibrium quantum field theory can be described by 
its classical field theory limit for sufficiently high initial 
occupation numbers. We show that eventually the cor- 
respondence breaks down due to the difference between 
classical and quantum thermal equilibrium. 

Employing the classical statistical field theory limit, 
we compare the 1/A^ expansion of the 2PI effective ac- 
tion at NLO with the exact result, which includes all 
orders in 1/A'^. The time evolution of classical nonequi- 
librium Green's functions can be calculated exactly, up 
to controlled numerical uncertainties, by integrating the 
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microscopic field equations of motion. This allows one 
to obtain a direct comparison ^. Apart from bench- 
marking approximation schemes employed in quantum 
field theory, the importance of the classical field limit for 
the approximate description of nonequilibrium quantum 
fields is manifest. 

2PI effective action. We consider a real A^-component 
scalar quantum field theory with a A((?!)a0a)^/(4!Af) in- 
teraction in the symmetric phase (a = 1,...,^^). The 
2PI generating functional for Green's functions can be 
parametrized as yj 

T[G] = ^TrlnG"^ + ^TrG^^G + T2[G] + const, (1) 

where Gq — i{0 + m?) denotes the free inverse prop- 
agator. The 2PI contribution r2[G] can be computed 
from a systematic 1/A^ expansion of the 2PI effective 
action [|. Writing T2[G] = T^°[G] + r^^°jp] -h ... the 
LO and NLO contributions are given by |^J^ 
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r^^"[G] = -Trcln[B(G)] 



(2) 
(3) 



Here C denotes the Schwinger-Keldysh contour along the 
real time axis iQ and 

B{x, y; G) - S^^+^x ^ v) + ^^ Gab{x, y)Gab{x, y). (4) 



The nonlocal four-point vertex at NLO is given by 
g^ B~^ p. In absence of extern 
equation for G is determined by 



g^ B ^ H . In absence of external sources the evolution 



6T[G] 
SGab{x,y) 



= 0. 



(5) 



We consider the 0{N) symmetric case, such that 
Gab{x, y) — G{x,y) Sab- In the following we solve Eq. (||), 
using Eqs. (M^M) without further approximations nu- 
merically in 1 -f 1 dimensions Q. We compare the out- 
come with results in the classical field theory limit, using 



both the NLO classical approximation and the "exact" 
Monte Carlo calculation ^ as described below. 

Nonequilibrium time evolution. To formulate the non- 
equilibrium dynamics as an initial- value problem, we de- 
compose the full two-point function using the identity 
G{x,y) = F{x,y) - (i/2)p(a;, y) signc(x° - y") where F 
is the symmetric or statistical two-point function and p 
denotes the spectral function [||. Following Refs. [§|,||, 
Eq. (^ can be written as 

p, + M^{x)] F{x,v)= - I dz"/"dz ^p{x,z)F{z,y) 

+ J dz°Jdzi:F{x,z)p{z,y), (6) 



dz°/ dz T,p{x,z)p{z,y). 



[U^ + M^{x)\p{x,y) 



At NLO in the 1/N expansion the effective mass term 
M'^{x) is given hYM'^{x) = m^ + X^^F{x,x) and the 
self energies are [|[ 

T.F{x,y)^-— F{x,y)lF{x,y) ~ -p{x,y)Ip{x,y) ,(7) 

Sp(a;,j/)=-— p{x,y)lFix,y) + F{x,y)Ip{x,y) . (8) 

Here the functions If and Ip resum an infinite chain of 
bubble diagrams, 

A A f^° f 

lF{x,y) = --IlF{x,y) + - / dz° dz I p{x,z)IlF{z,y) 

A f^" f 

A A f''^° f 

Ip{x,y) == --Tlp{x,y) + - / dz° / dz Ip{x,z)'^p{z,y), 

with 

^f{x. y) = -\{f^{^-^ v) - -/{^^ y)) : (9) 

np{x,y)^-Fix,y)p{x,y). (10) 

Classical field theory limit. The classical statistical 
field theory limit of a scalar quantum field theory has 
been studied extensively in the literature. An analysis 
along the lines of Refs. ||ll|-[l4[| shows that all equations 
(m)-(10[) remain the same in the classical limit except for 
differing expressions for the statistical components of the 
self energy 
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(11) 
(12) 



The latter expressions are given by Eqs. (0) and (^ in 
the quantum theory. One observes that the classical self 



energies are obtained from the expressions in the quan- 
tum theory by dropping terms with two spectral (p-type) 
components compared to two statistical (F-type) func- 
tions. This particular relationship has been studied in 
great detail in thermal equilibrium and corresponds to 
retaining only contributions that are of leading order in 
h. It has been systematized in terms of Feynman rules 
for classical and quantum theories using the Keldysh for- 
mulation with appropriate interaction vertices [Q . Note 
that the classical spectral function is obtained from the 
quantum one by replacing —i times the commutator by 
the classical Poisson bracket. A comparison of the clas- 
sical limit in the current approximation for TV = 1 has 
been studied in Ref. Q. 

Monte Carlo approach. An "exact" nonperturbative 
solution of the evolution of classical correlation functions 
in the 0{N) model can be obtained numerically in a 
straightforward manner B . Initial conditions are deter- 
mined from a probability functional on classical phase- 
space. The subsequent time evolution is solved numer- 
ically using the classical equations of motion. In the 
figures presented below, we have sampled 50000-80000 
independent initial conditions to approximate the exact 
evolution. 

Far from equilibrium evolution. We consider a system 
that is invariant under space translations and work in 
momentum space. We choose a Gaussian initial state 
such that a specification of the initial two-point func- 
tions is sufficient. The quantum (classical) spectral 
function at initial time is completely determined from 
the equal-time commutation relations (Poisson brack- 
ets). For the symmetric two-point function we take 
F{0,0;p) = [no{p) + ^]/u!p, with the initial particle 
number no(p) — ritsip) + nsip) representing a peaked 
"tsunami" nts{p) — Aexp [~2^i\p\ ~ bts|)^] in a ther- 
mal background 7is(p) = [exp(a;p/To) — 1]^^ ||l5|,p|,p| . 
Such an initial state is reminiscent of two colliding wave 
packets |15|| and provides a far-from-equilibrium initial 
condition. We emphasize that these initial conditions can 
be implemented both in a quantum and a classical theory. 
The initial mode energy is given by ujp = {p^ + Af^)^/^ 
where M is the one-loop renormalized mass in pres- 
ence of the nonequilibrium medium, determined from the 
one-loop gap equation. As a renormalization condition 
we choose the one-loop renormalized mass in vacuum 
TOfl = M|„(,=o = 1 as our dimensionful scale. The re- 
sults shown below are obtained using a fixed coupling 
constant A/m|j — 30. 

Convergence of NLO and Monte Carlo results. In 
Fig. n we present the two-point function F{t,0;p = 0) 
in the classical field theory limit for three values of N. 
All other parameters are kept constant. The figure com- 
pares the time evolution using the 1/iV expansion of the 
2PI effective action to NLO and the Monte Carlo cal- 
culation that includes all orders in 1/A^. One observes 
that the approximate time evolution of the correlation 
function shows a rather good agreement with the exact 
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FIG. 1. Unequal time two-point function F{t,0;p — 0) at 
zero momentum times the initial mass M for N — 2, 10, 20. 
The full lines show results from the NLO classical evolution 
and the dashed lines from the exact classical evolution (MC). 
For A'^ = 20 the NLO and exact evolution can hardly be dis- 
tinguished. 

result even for small values of N (note that the effective 
four-point coupling is strong, X/6N — 2.5to|j, for N — 2). 
For A^ = 20 the exact and NLO evolution can hardly be 
distinguished. A very sensitive quantity to compare is 
the damping rate 7, which is obtained from an exponen- 
tial fit to the envelope of F{t,0;p = 0). The systematic 
convergence of the NLO and the Monte Carlo result as a 
function of 1/A^ can be observed in Fig. g. The quanti- 
tatively accurate description of far from equilibrium pro- 
cesses within the NLO approximation of the 2PI effective 
action is manifest. 

Classical behavior of nonequilibrium quantum fields. In 
Fig. 1^ we also show the damping rate from the quantum 
evolution, using the same initial conditions and parame- 
ters. We observe that the damping in the quantum the- 
ory differs and, in particular, is reduced compared to the 
classical result. In the limit N ^ 00 damping of the 
unequal-time correlation function F(t,0;p) goes to zero 
since the nonlocal part of the self energies (0)-(P van- 
ishes and scattering is absent. In this limit there is no 
difference between evolution in a quantum and classical 
statistical field theory. 

For finite TV scattering is present and quantum and 
classical evolution differ in general. However, the clas- 
sical field approximation may be expected to become a 
reliable description for the quantum theory if the num- 
ber of field quanta in each field mode is sufficiently high. 
We observe that increasing the initial particle number 
density leads to a convergence of quantum and classical 
time evolution at not too late times. In Fig. 0we present 
the time evolution of the equal-time correlation function 
F{t,t;p) for several momenta p and N — 10. Here the 



FIG. 2. Nonequilibrium damping rates extracted from 
F{t,0;p — 0) shown in Fig. hi as a function of 1/A''. Open 
symbols represent exact and NLO classical evolution. One 
observes a rapid convergence of the 1/N expansion at NLO 
to the exact MC result. The quantum results are shown with 
full symbols. In the quantum theory the damping rate is re- 
duced compared to the classical theory. 

particle density / ■^nQ{p)/Ad — 1.2 is six times as high 
as in Figs. wM and, in contrast to the latter, quantum 
and classical evolution at NLO follow each other rather 
closely. For an estimate of the NLO truncation error we 
also give the MC result for A^ = 10 showing a quantita- 
tive agreement with the classical NLO evolution both at 
early and later times. 

Quantum versus classical equilibration. From Fig. y 
one observes that the initially highly occupied "tsunami" 
modes {pts/'mR = 2.5) "decay" as time proceeds and 
the low momentum modes become more and more pop- 
ulated. At late times the classical theory |1^J^ and the 
quantum theory ^,g] approach their respective thermal 
equilibrium distribution. Since classical and quantum 
thermal equilibrium are distinct the classical and quan- 
tum time evolutions have to deviate at sufficiently late 
times, irrespective of the initial particle number density 
per mode. Differences in the particle number distribution 
can be conveniently discussed using the inverse slope pa- 
rameter T{t,p) = —n{t,ep)[n(t,ep) + l]{dn/de)~^ for a 
given time-evolving particle number distribution n{t, e, 
and dispersion relation ep{t) g. Following Ref. 
we define the effective particle number as n{t, e^) + 
\ = [F{t,t';p)dtdt'F{t,t';p)Y/\=t' and mode energy 
by ep{t) = [dtdt'F{t,t'-p)/F(t,t'-p)Y'^\t=t', which coin- 
cide with the usual free-field definition for A ^ 0. For 
a Bose-Einstcin distributed particle number the parame- 
ter T{t,p) corresponds to the (momentum independent) 
temperature T{t,p) = Teq. In the classical limit the in- 
verse slope T{t,p) as defined above remains momentum 
dependent. 
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FIG. 3. Nonequilibrium evolution of the equal-time 
two-point function F{t,t;p) for A'^ = 10 for various momenta 
p. One observes a good agreement between the exact MC 
(dashed) and the NLO classical result (full). The quantum 
evolution is shown with dotted lines. The initial particle den- 
sity is six times as high as in Figs. hlM At these high densi- 
ties, the difference between quantum and classical evolution 
is small. 

In Fig. ^ we plot the function T{t,p) for piow — and 
Phigh — 2pts- Initially one observes a very different be- 
havior of T(t,p) for the low and high momentum modes, 
indicating that the system is far from equilibrium. Note 
that classical and quantum evolution agree very well for 
sufficiently high initial particle number density. However, 
at later times the difference between quantum and clas- 
sical evolution becomes visible. The quantum evolution 
approaches quantum thermal equilibrium with a momen- 
tum independent inverse slope T — A.lmfi. In contrast, 
in the classical limit the slope parameter remains momen- 
tum dependent and the system relaxes towards classical 
thermal equilibrium |0|. 
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